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We report the first observation of i?^ — > '>p{2S)<l) decay in pp collisions at 
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^/s = 1.96 TeV using 360 pb of data collected by the CDF II detector at the 
Fermilab Tevatron. We observe 20.2 ± 5.0 and 12.3 ± 4.1 S° ^ V'(25')0 candidates, 
in ip{2S) l-i^fJ'" and '4'{2S) J / iIjti'^ -n" decay modes, respectively. We present a 
measurement of the relative branching fraction B{B^ ilj{2S)(j))/B{B'^ — > J/V"/*) = 
0.52 ± 0.13(stat.) ± 0.04(syst.) ± 0.06(BR) using the ^'(25') fi~ decay mode. 

PACS numbers: 13.25.Hw 



The decays of B mesons to charmonium final states have been studied extensively in the 
past, and the measurements Q, 0] show that the ratio of the branching fractions of 
and decay to the ip{2S) final states over the J/ip final states are approximately 60% 
as shown in Table E The 5+'° ^ J/ipK+'*^ {ip{2S)K+'*°) and 5° J/ifjcf) (^(2^)0) are 
color-suppressed Cabibbo-favored decays that have the same tree-level decay topology as 
shown in Fig^ The relative branching ratio between B^ — > ■il){2S)(f) and 5° J/i^P has 
not been measured. Only one 5^ — ^ ilj{2S)(f) candidate event has been reported at LEP in 
1993 Q. 



b;,b°,b° 



c 

J/V, \|/(2s) 

c 



u, d, s u, d, s 



FIG. 1: Tree level Feynman diagram of B mesons decaying to charmonium final states. 



TABLE I: The current relative branching ratio of B meson decays between ip{2S) and J ftp final 
states. 



Decay channel Value Reference 



B{B+^J/^K+) 0.64 ± 0.06 ± 0.07 BaBar QJ 



B{B+^i;{2S)K+) 
B+^J/^ipK+) 

B(B^^ij(2S)K^) 

B{B^Zj/^K^) 0.82 ± 0.13 ± 0.12 PDG U 
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The J/'ip4> mode has recently been used to determine the decay widths for the 

heavy and hght mass eigenstates by measuring the relative contribution of the CP-odd 
and CP-even components to the observed angular distribution as a function of the decay 
time 0, Q|- Observing the B^ il){2S)(j) would allow an independent measurement of the 
decay widths for the heavy and light mass eigenstates in the future. In particular, the 
spin alignment of 5° il){2S)(j) could be different from that of 5^ J/tpcp. 

In this Letter, we report the observation of B^ il){2S)(j) in both ip{2S) fi'^fi~ and 
^^{23) J/ipTx^Ti" decay modes produced in pp collisions at ^/s = 1.96 TeV. We also mea- 
sure the ratio of branching fractions for B^ J/'ip4' B^s ~^ i'{'^S)(f). Many systematic 
effects cancel in the measurement of the ratio, including uncertainties in total integrated 
luminosity, bottom-quark production and fragmentation, and trigger and reconstruction ef- 
ficiencies. In addition, for this ratio of branching fractions measurement we use only the 
'0(25') — > fJ'~^fJ'~ decay mode in order to guarantee identical topologies for the J/tp and ip{2S) 
channels. Similar decay modes, such as B^ J/tpK^ and B^ tlj{2S)K~^, are used as con- 
trol samples to perform consistency checks, and study the systematic uncertainties. Charge 
conjugate modes are implied throughout this paper. The data sample is comprised of about 
3 X 10^ J/ij fi+fi-, 1 X 10^ ^ip{2S) ii+fi-, and 1.6 x 10^ ij{2S) J/^Ti+n- candi- 
dates. The total integrated luminosity is approximately 360 pb~^ and was collected using 
the Collider Detector at Fermilab (CDF II) between February 2002 and July 2004. 

The CDF II detector is described in detail elsewhere P|. The main components for this 
analysis are tracking and muon systems. The tracks are reconstructed by the silicon mi- 
crostrip detector (SVX II) j^, the intermediate silicon layer (ISL) [l^, and the Central Outer 
Tracker (COT) [111, which are immersed in a uniform axial 1.4 T magnetic field provided 
by a superconducting solenoid. The SVX II is located at a radial distance between 2.5 
and 10.6 cm, with double-sided micro-strip sensors arranged in five concentric cylindrical 
shells. The ISL provides one additional layer for central tracks, and two for tracks with 
|?7| > 1.0. The layers are located radially between the SVX II and the COT. The COT is 
a multi-wire drift chamber, which consists of 96 layers of sensitive wires, grouped into 8 
superlayers. These superlayers cover the radial region between 40 cm and 137 cm from the 
beam, and are arranged to provide alternating axial and ±2° tilted stereo measurements. 
Planar drift chambers ^ | located outside the calorimeter are used to identify muons in the 
central region (|?7| < 1.0, where r] is the pseudorapidity) . The events are selected with a 
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three-level trigger system. At Level 1, charged particle trajectories in the plane transverse 
to the beam direction are reconstructed from the COT hits using a hardware processor jisl ]. 
The trigger requires tracks with transverse momentum Pt(/^) > 1-5 GeV/c to be matched to 
hits in the muon detector. At Level 2, opening angle and opposite-charge cuts are imposed 
on the muon pairs. At Level 3, the two muon tracks are required to be oppositely charged 
with invariant mass between 2.7 and 4.0 GeV/c^. 

We reconstruct 5° ^ and 5° ^ ^(2^)0 followed by ip{2S) fx+fj.- and ip{2S) 

J/ipTT^TT^, where J/ip ^ fi'^fi" and — s> K^K~. For the measurement of the relative 
branching fraction between i?^ J/i^P and — > ■ip{2S)(j), it is desirable to have selection 
criteria similar for both decay modes. All three decay channels involve only the well 
known J/ip, '0(25') and decays, which have been used extensively in other measurements 
at CDF, and their selection criteria are well established. In this analysis, we follow the 



14| and apply them 



selection requirements developed in the b hadron mass measurements 
to the three decay modes of interest. 

The reconstruction begins by selecting J /ip ^ or '0(25') fi~^fi~ candidates, with 

pairs of oppositely charged tracks that satisfy the muon pair trigger requirements. The 
reconstructed /i^yU" invariant mass is required to be within 80 MeV/c^ of the J ftp or 0(25) 
mass P]. The ^^{28) — J/ipn'^'n" is reconstructed by associating a. J/ip ^ i^i'^fi~ candidate 
(with its mass constrained to the J/ip mass) with a pair of tracks, each with > 0.4 GeV/c. 
The invariant mass of J/ipir^ir" is required to be within 20 MeV/c^ of the world average 
■0(25) mass ^. Once a J/tp or 0(25) candidate is selected, we search for a — > K^K~ 
candidate with a pair of additional tracks. The invariant mass of K^K~ is required to be 
within 10 MeV/c^ of the mass p|. The pt of the candidate is required to be greater 
than 2.0 GeV/c. The meson candidates are then reconstructed by associating a J/tp or 
'0(25') candidate with a candidate. All tracks (4 tracks in B^ J/ip<P or 5° '0(25')0 
followed by 0(25) and 6 tracks in i?^ ■ip{2S)(j) followed by 0(25) J/ipn^TX") 

are required to be consistent with having originated from a common vertex satisfying vertex 
quality requirements. Prompt background, with tracks coming directly from the primary 
vertex, can be reduced by exploiting variables sensitive to the long lifetime of the meson. 
To reduce prompt background, the transverse decay length {L^y) of the 5° is required to 
exceed 100 /im, where L^y is defined as the transverse vector from the beam axis to the 
B^ decay vertex projected onto the transverse momentum of the 5° candidate. To ensure 
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a well measured B meson decay vertex, each track is required to have a measurement in 
at least three axial layers of the silicon detector, including SVXII and ISL. The transverse 
momentum of the candidates is required to be greater than 6.5 GeV/c to further reduce 
combinatoric background. To improve the B meson mass resolution, the fi'^fj^" mass is 
constrained to the J/ip or ip{2S) mass, while the J /ipir^ii' mass is constrained to the ip{2S) 
mass. 
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5.3 5.35 
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FIG. 2: Invariant mass distributions for — > J/ipcj) (bottom), and for B'^ 'tp{2S)(l), followed 
by ip{2S) — > fi^ ^~ (upper left), or il^{2S) J/V'vr+vr" (upper right). The curves are the results 
of the fits described in the text. 



Two sources of background are expected in the 5° signal region: combinatoric background 
and kinematic "reflection" of 5° ^ J/i)K*^ (for 5° ^ J/^0) or 5° ^ i){2S)K*^ (for 5° ^ 
tlj{2S)(f)), where the pion from the K*^ decay is mis-assigned as a kaon. The combinatoric 
background is modeled by a first order polynomial. The B^ — * J/ipK*^ {B^ il){2S)K*^) 
reflection background results in a broad distribution near and above the signal region. 
The fraction of 5° ^ J/^jK*^ {B° ^{2S)K*^) events that fah into the 5° ^ J /ijcj) 
(5° — s> il){2S)(j)) signal region is estimated using a Monte Carlo simulation. The background 
contribution from reflection in our data sample is then calculated by multiplying the fraction 
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determined from Monte Carlo simulation by the number of the J/ipK*^ and 

i){2S)K*^ candidates in the same data. The contribution of the B^ J / ipK*^ reflection in 
the 5° J/'4'4' signal region is estimated to be 6.6 ± 0.3 events. The contribution of the 
B^ il){2S)K*^ reflection in the i?^ ■il){2S)(f) signal region is estimated to be 0.34 ± 0.05 
and 0.19 ± 0.03 events for ^/'(2S') yU"*"/^" and ip{2S) J/tp7r'^7r~ modes, respectively. 
The B^ J/tpK*^ [B^ il){2S)K*^) reflection background is highly suppressed because 
only a small fraction of the misidentified K*^ — > A'+vr" can satisfy the cf) K^K~ mass 
requirement. 

TABLE II: The numbers of observed signal events, the fitted masses, and the signal width (Gaussian 
sigma, or mass resolution) for each of the three decay channels. Note that the width for each 
of ip{2S)4> decay modes is fixed as described in the text. 



Decay Mean[MeV/c2] Width[MeV/c2] Yield 

5° ^ J/i^cj) 5366.76±0.66 9.42±0.58 292.2± 15.9 

5° ^ 'i/'(25')(/>;'0(25') ^ 5366.50±1.86 6.63 (fixed) 20.2 ±5.0 

5° ^ ^(25)(/>;'0(2S') ^ J/V'vr+Tr- 5366.63±3.20 7.77 (fixed) 12.3±4.1 



An unbinned log-likelihood fit is used to extract signal yields from the reconstructed 
mass spectra, as shown in Fig. |21 The signal distribution is modeled as a Gaussian, and the 
background distribution is modeled as a first order polynomial. The background component 
from misidentified K*^ decays is also included, with a shape obtained from the Monte Carlo 
simulation. The width for each of the two 5° il){2S)(j) modes is fixed in the following 
way. We take the ratio of the widths for i?^ il){2S)(j) relative to i?^ J/i"P as a 
scale factor determined from Monte Carlo simulation, and we then calculate the width for 
-B^ — > tlj{2S)(f), using the width of B^ J/i"P from data. Comparison between Monte 
Carlo and data for the control samples of B^ J/ipK'^ and B^ il){2S)K^ shows that 
the relative ratio of the widths of the two modes can be well predicted by Monte Carlo 
simulation. The signal yields, fitted masses, and width (Gaussian sigma, or mass resolution) 
of the three decay channels are summarized in Table |nj A consistency check (Monte Carlo 
independent) is performed by fitting the (f) K^K~ invariant mass spectra for events in 
the Bs signal region after sideband subtraction. The K'^K~ signal yield obtained this 
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way in data is consistent with that from fitting the Bs mass spectra, indicating that the 
reflection background from B^ J/tpK*^ {B^ — > ip{2S)K*^) or other decay modes is 
indeed neghgible.. 

The background contribution in the signal region (defined as a window six times the 
expected mass resolution, as shown in Table II, around the mean value of the Bs signal 
peak) for 5° ip{2S)(j), followed by ipi^S) fJ'~^fJ'~ and ip{2S) J/ip7r^7r~ decays, is 
estimated to be 10.0 ± 3.2 and 6.5 ± 2.6 events, respectively. The probability of a statistical 
fluctuation of the expected total background in the signal region to the observed or higher 
number of events is 2.5 x 10"'' for B^ il){2S)(j) with ip{2S) —>■ fi^fj^~ and 1.6 x 10~^ for 
ip{2S) — * J/ipTT^ir". These correspond to 5.0 cr and 4.2 cr one-sided Gaussian significance 
for the two decay modes, respectively. The combined probability of the two modes is 1.1 x 
10"^'^, corresponding to a 6.4 a significance for the observation of B^ il){2S)(f). 

We measure the relative branching fraction between 5^ J/'4^'P and -B^ tp{2S)(f) 
using only the J/tp ^ f^~^fJ'~ decay mode, and the control sample data (5+ J / tpK^ and 
B^ iIj{2S)K^) are used to study the systematic uncertainties. The relative branching 
ratio for the fi~^fJ^~ mode is extracted using the formula: 



BjB^s ^ ^(2^)0) _ iV^(25)^ BjJ/jj ^ ej/^^ 

where ej/^^/e^(2S)<i> = 0.925±0.006 is the ratio of the combined trigger and selection efficien- 
cies derived from Monte Carlo simulation (with the error due to the size of the simulated 
samples), and Njj^^ or N^(2S)^ is the total number of reconstructed B^ mesons for each 
mode. The B{J/ip — /i'^/i~) and B{ip{2S) f^~^f^~) are the world average branching 
fractions 0. 

In our analysis, we use a Monte Carlo simulation to determine the relative efficiency 
for the two decay modes, and the control sample data are used to study the systematic 
uncertainties. The simulation of the CDF II detector is based upon a geant description jlS^ . 
Transverse momentum and rapidity distributions of single b quarks are generated based on 
next-to-leading-order (NLO) pertubative QCD The B^ meson spectrum used in the 

Monte Carlo simulation is consistent with the data from inclusive B J/ipX jsl. The 
EVTGEN program ^ is used to decay B mesons into the final states of interest. 

Since both modes are B^ decays, and the decay topologies are very similar, most sys- 
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tematic uncertainties cancel in the ratio. Systematic uncertainties originate from fitting the 
invariant mass distributions to obtain signal yields, from determination of the relative effi- 
ciencies, and from the measured branching fractions of J/ip and ipi'iS) decays taken from j^. 
Consistency checks are performed on the fitting method by varying the range and using dif- 
ferent functions, and no statistically significant variation is found. Systematic uncertainty 
from the fitting method is evaluated by performing the fitting without constraints on the 
width. 

Systematic uncertainties on the ratio of efficiencies are due to the differences in the 
kinematics of the two decay modes. For example, due to the mass difference between ip{2S) 
and J/ip, the prif^) distributions are somewhat different between the two decay modes. To 
take into account the difference in prifJ') distributions, the single muon efficiency measured 
from data jQ] is used to reweight the Monte Carlo samples, and the relative efficiency 
(central value and error) is recalculated. We vary the measured muon efficiency, and find 
that the systematic uncertainty due to the difference in prifJ') distributions is negligible. The 
main systematic uncertainty due to decay kinematics difference comes from lack of knowledge 
of the angular correlation in the il){2S)(f) decay. The central value of the relative 

efficiency is determined by assuming that the angular correlation of the -B^ — > ■il){2S)(f) decay 
is the same as that of the 5^ J/ipcj). To evaluate the effects on our measurement, we 
generate Monte Carlo samples with pure CP-even and CP-odd decays for i?^ ?/;(25')0 
and recalculate the relative efficiency. We take the difference between CP-even and CP-odd 
cases as the systematic uncertainty, which turns out to be the major component (5.5%). The 
systematic uncertainty from the fitting contributes at the 3.9% level. The total systematic 
uncertainty is 6.7%. 

The contribution from the branching fractions is calculated by propagating the world 
average uncertainties. The dominant contribution is due to the measured branching ratio 
B{ij{2S) ^ ^+/i-) = (0.73 ± 0.08)% 0. 

Using Eq. ^ we derive the ratio of relative branching fractions: 

'^fnt^ = 0-52 ± 0.13(.tal) ± QM{syst.) ± 0.06(5/?) (2) 
where the first error is statistical, the second is systematic, and the third is due to the 
branching ratios oi J/ip f^~^f^~ and ip{2S) /^"''yU"- 

In summary, we present the first observation of 5° ip{2S)(j) decay, in both iIj{2S) — >• 
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/i"*"/!" and ip{2S) J/ipix^Ti" modes mpp collisions at y/s = 1.96 TeV using the CDF II de- 
tector. We also present the measurement of the ratio of branching fractions between 
Bg '^(25')0 and 5° J/4''P using the '1^(25') fi~^fi~ decay mode. This result for 

is consistent with the ratios of branching fractions for the corresponding decays of 
and B^ P], indicating that the relative branching ratio of B meson decays between ip{2S) 
and J/ ip final states is independent of the flavor of the lighter quark. 
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